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Abstract 

Recent measurements of PeV energy neutrinos at IceCube and a 3.5 keV X-ray line in the spectra 
of several galaxies are both tantalizing signatures of new physics. This paper shows that one or both 
of these observations can be explained within an extended supersymmetric neutrino sector. Obtaining 
light active neutrino masses as well as phenomenologically interesting (keV-GeV) sterile neutrino 
masses without any unnaturally small parameters hints at a new symmetry in the neutrino sector 
that is broken at the PeV scale, presumably tied to supersymmetry breaking. The same symmetry 
and structure can sufficiently stabilize an additional PeV particle, produce its abundance through the 
freeze-in mechanism, and lead to decays that can give the energetic neutrinos observed by IceCube. 
The lightest sterile neutrino, if at 7 keV, is a non-resonantly produced fraction of dark matter, and can 
account for the 3.5 keV X-ray line. The two signals could therefore be the first probes of an extended 
supersymmetric neutrino sector. 
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I. INTRODUCTION AND MOTIVATION 


Astrophysical signatures of dark matter (DM) interactions, coupled with theoretical and 
observational input from other areas of particle physics, can provide crucial insights into the 
nature of physics beyond the Standard Model (BSM). However, while the existence and macro¬ 
scopic properties of dark matter in our Universe have been well established, its microscopic 
properties continue to remain elusive. Compelling theoretical arguments dictate that dark 
matter is most likely composed of weakly interacting massive particles (WIMPs), predicting a 
GeV-TeV mass scale for dark matter, although supporting observational evidence has so far 
failed to materialize. Interestingly, some recent observations hint at dark matter beyond this 
narrow GeV-TeV window, indicating that the structure of dark matter and BSM physics might 
be very different from what was envisioned. 

One such observation is the IceGube neutrino observatory’s detection of 37 neutrino events 
in the energy range from 30 GeV to 2 PeV, disfavoring a purely atmospheric explanation at 5.7 ct 
PP. In particular, IceGube has reported 3 neutrino events above 1 PeV. Strongly incompatible 
with traditional astrophysical processes, these have been shown to be compatible with the decay 
of heavy, long lived particles with lifetimes of about lO^^s that constitute some or all of dark 
matter [2HI2]. A dark matter interpretation of these events, however, poses several theoretical 
inconveniences: Why should we expect new physics at the PeV scale? Why should a PeV scale 
particle be so long-lived? What sets its relic density? Why should it show up in neutrinos? 
Given the dearth of BSM candidates at the PeV scale, a theoretically motivated framework 
that also hts into the broader particle physics picture would be extremely appealing. 

Another observation, made independently by two groups, is the discovery of an unidentified 
X-ray line at ^ 3.5 keV in the stacked X-ray spectra of 73 galaxy clusters measured by 
XMM-Newton PI and in the X-ray spectra of the Andromeda galaxy and the Perseus galaxy 
cluster HI. A dark matter interpretation of these observations is plagued by many unresolved 
questions, such as possible contamination from potassium lines HMD as well as inconsistency 
with stacked observations of dwarf spheroidals [18] and galaxy spectra [19] (although greater 
compatibility is reported with observations of the Milky Way [201^1]). Despite these concerns, 
we entertain the possibility that these observations can be explained by a monochromatic line 
from the decay of a 7 keV dark matter particle. Unlike the IceGube PeV neutrinos, there exists 


2 


an extremely well-motivated keV scale dark matter candidate in the form of a sterile nentrino, 
whose decay into an active nentrino and an X-ray photon with energy eqnal to half its mass 
has long been heralded as its smoking gnn signatnre. Sterile nentrinos are essential ingredients 
of seesaw models of nentrino masses, and while their masses conld, a-priori from theory, lie 
anywhere from the eV to the GUT scale, forming a significant fraction of dark matter while 
remaining consistent with X-ray [22H2S] and Lyman-alpha [HOHM] measnrements constrains 
them to be at the keV scale. 

The pnrpose of this paper is to show that one or both of these signals can naturally arise 
from an extended supersymmetric neutrino sector. The framework is a simple extension of a 
recently proposed model [3Z] , where it was shown that if the right-handed neutrinos are charged 
under some new symmetry U(l)', broken by the PeV scale vacuum expectation value (vev) of 
a scalar field, active neutrino masses consistent with data and keV-GeV mass sterile neutrinos 
that can form part or all of dark matter can be realized without the need for any unnaturally 
small parameters in the theory. The PeV scale vev can presumably be inspired by the scale of 
supersymmetry breaking. There exist compelling arguments for supersymmetry at such high 
scales from flavor, GP, and unification considerations High scale SUSY is now further 

motivated by null results at numerous indirect and direct detection experiments and the first 
run of the LHG. The Higgs boson mass rrih = 125 GeV is also compatible with PeV scale 
sfermions for 0(1) values of tan/? in the MSSM 

The realization of the 3.5 keV line from a sterile neutrino component of dark matter in this 
framework is straightforward. In addition, we will see that a straightforward extension of the 
neutrino sector using the same symmetry and structure employed for neutrino masses 

can result in a PeV scale dark matter candidate whose decays are compatible with the high 
energy neutrino events at IceGube. It is interesting to note that a mixture of cold and warm 
dark matter components might also be favorable for solving some of the small scale problems 
in cosmology [IHl Il6] . 

Finally, a few words on the philosophy behind the structure paper. This paper does not 
aim to build a complete model of PeV scale supersymmetry - this would bring in model- 
dependent details and complexities that are irrelevant to the task at hand. We simply intend 
to demonstrate that our framework can naturally accommodate the IceGube PeV neutrinos and 
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the 3.5 keV X-ray line; for this reason, we also do not pursue detailed scans over parameters, 
and the benchmark values presented in this paper should not be interpreted as the “best-£t” 
ones. The greatest virtue of this exercise lies not in obtaining best fits or building a complete 
model but in placing the two signals in perspective within a broader particle physics framework, 
where they are not isolated observational curiosities but are connected to outstanding issues in 
the neutrino and Higgs sectors. Finally, although unifying such disparate scales as the keV and 
PeV into a single simple framework appears extremely appealing, it should be kept in mind 
that these two signals are by no means necessary ingredients in the theory; indeed, both, one, or 
neither of them can be realized in the model with appropriate choices of helds and parameters. 


II. THE MODEL 

In this section we present our model of the supersymmetric neutrino sector, which aims to 
realize active neutrino masses consistent with oscillation data, a keV sterile neutrino dark matter 
candidate, and a PeV scale dark matter candidate without any unnaturally small parameters 
in the theory. This is an extension of the model in [37j, and the interested reader is referred 
to this paper for further details. As motivated in B7|, we take the scale of supersymmetry 
breaking to be around 0{1 — 100) PeV. This means that all dimensionful parameters obtained 
after supersymmetry breaking, such as masses and vevs, are expected to be at this scale. 

In order to obtain neutrino masses, three SM-singlet sterile neutrinos Ni are introduced. 
Although the Ni are singlets under the SM gauge group, they are unlikely to be singlets under 
all symmetries of nature, as this would naturally place their masses at the Planck scale or 
the GUT scale, contrary to what is phenomenologically desirable. We therefore posit that the 
Ni are charged under some new symmetry of nature (for concreteness, a U{!)'). The W are 
coupled to a single exotic field 0 of equal and opposite U(l)' charge to form singlets, 

which can then be coupled to SM fields. 

With the IceCube PeV neutrinos in mind, we also introduce a new held X, a SM singlet 
scalar with a PeV scale mass and appropriately charged under the U(l)' to be sufficiently long- 
lived to form a component of dark matter. It is likewise coupled to another held Y (which 
carries lepton number) to form a U(l)' singlet. Given the new symmetry U(l)', this is the most 
straightforward extension of the model to incorporate additional helds. 
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Supermultiplet 

spin 0, 1/2 

u{iy 

Remarks 

Mi 

Ni, Ni 

+1 

Ni sterile neutrinos 


4>,y’4‘ 

-1 

((/>) ~PeV, breaks 17(1)' 

X 

X, V’x 

+5 

mx ~PeV, dark matter 

y 

Y, 

-5 

U{iy partner of X 


TABLE I. Field content, notation, and U{\y charge assignments for the new multiplets introduced in 
the neutrino sector of the model. These lead to the higher-dimensional operators in the superpotential 
in Eq.[TJ 

Since the theory is supersymmetric, each of these helds resides in a chiral supermultiplet; 
the held content and notation are summarized in Table These lead to the following non- 
renormalizable terms in the superpotential that are relevant to our study: 


M* M* M* M* M* 


1 -^2 , 1 ^3 


5! Mi* 


+ 


5! M3 


yyfl 


( 1 ) 


All couplings are written as dimensionless numbers and expected to be 0{1). The iVj basis is 
chosen such that the third term in Eq.[^ is diagonal. M* is the scale at which this effective 
theory of non-renormalizable operators needs to be UV completed with new physics, such as 
the scale of grand unihcation Mgut or the Planck scale Mp. 

Obtaining Dirac and Majorana masses for the sterile neutrinos Ni in order to recover the 
seesaw mechanism requires the scalar component 0 of <h to obtain a vev, thereby breaking the 
U{1)' symmetry. We assume that (j) obtains a PeV scale vev from the supersymmetry breaking 
sector (without delving into details of how exactly this might be realized, which is tangential 
to the main purpose of this paper). In addition, we also assume that the helds in the X,y 
and $ multiplets all get PeV scale masses. This setup has the following phenomenological 
consequences: 
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A. Neutrino Masses 


With <p obtaining a vev at the PeV scale and Hu acqniring a vev from electroweak symmetry 
breaking, the hrst and third terms in the superpotential in Eq.[T]lead to the following active- 
sterile Dirac masses and sterile Majorana masses in the neutrino sector (flavor indices suppressed 
for simplicity): 


C(0)(g°) 


rriM 


M, 


( 2 ) 


The seesaw mechanism then gives the following sterile and active neutrino mass scales: 

ml CHHl? 


ms = tum = 


M* 


rua = 


rriM 


rfM^ 


( 3 ) 


which also determines the mixing angle between the active and sterile sectors: 

cm 


9^ J— = 


m^ 




( 4 ) 


With (0) ~ 1 -100 PeV, M, = Mgut = 10^® GeV, tan/3 = 2 {{H^) = 155.6 GeV), and 0{1) 
values of ( and 77 , this framework produces active neutrino masses that £t oscillation data 
and sterile neutrinos with (P(keV-GeV) masses, which are compatible with dark matter and 
cosmological observations (see 133 for more details). 


B. Sterile Neutrinos and Dark Matter 

The three sterile neutrinos W naturally have masses at the keV-GeV scale in this framework. 
We require the lightest one, Ni, to be a dark matter candidate with a keV scale mass in order 
to explain the 3.5 keV X-ray signal. However, several recombination era observables 
constrain the two heavier sterile neutrinos ^ 2,^3 to decay before Big Bang Nucleosynthesis 
(BBN), forcing Tm,N 3 ^ Is and consequently mm ,ns ^ (P(IOO) MeV. This mass hierarchy 
between m^i and mm ,3 requires a similar hierarchy between (ij and rji values, necessitating 
some tuning of parameters. 

As Ni couples extremely weakly to the SM fields and is never in thermal equilibrium in 
the early Universe, its relic abundance is not set by thermal freeze-out. It is produced instead 
through active-sterile oscillation at low temperatures, known as the Dodelson-Widrow (DW) 
mechanism [50]. A combination of X-ray bounds and Lyman-alpha forest data now rule out 
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the prospect of all of dark matter consisting of Ni produced via the DW mechanism; however, 
it can still constitute a signihcant fraction of the dark matter abundance (this will be further 


discussed in Section IIIB). This is desirable for us, since the remaining dark matter can come 
from the PeV sector. 


C. PeV Scale Dark Matter 

Our PeV dark matter candidate is a SM singlet scalar X that carries a charge (see 

Tablej^. In general, a PeV scale dark matter candidate presents two caveats. First, obtaining 
the correct relic density through the well-known thermal freeze-out mechanism requires a large 
annihilation cross-section that comes into conflict with unitarity limits [51], so that freeze out 
of a PeV scale particle overdoses the Universe. Hence one has to ensure that the candidate 
is never in thermal equilibrium and build up its abundance through some other mechanism. 
Second, since the decay rate of a particle is in general proportional to its mass, an unstable PeV 
scale particle is generically far too short-lived to be a dark matter candidate, and appropriate 
measures need to be put in place to stabilize it against rapid decay. We will see that the new 
U{iy symmetry in our theory can be used to address both of these issues. 

In light of the unitarity bound mentioned above, some further assumptions need to be made 
about the superpartners that freeze out of the thermal bath, as we take supersymmetry to be at 
the PeV scale. One possibility is to make the lightest supersymmetric particle (LSP) sufficiently 
light that the thermal freeze-out abundance does not overdose the Universe. Another is to have 
all superpartners decay into SM states through R-parity violating (RPV) interactions. In this 
paper we choose the former option despite the dear need to tune parameters in order to achieve 
this, as the latter would involve the introduction of several RPV operators in our model, which 
require careful treatment beyond the scope of this work. To this end, the LSP is chosen to be 
a Higgsino at ~ 800 GeV, which would make up about half of the dark matter abundance. 

1. Production of X 

Since X is charged under the U(l)', there are no renormalizable terms in the superpotential 
that connect it with SM fields. The lowest dimension term allowed, which must be both a 
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SM singlet and a U{iy singlet, is the term -^LiH^Xy (see Eq.j^; this leads to the following 
production processes for X from the thermal bath: 

Ih^Xi^Y, IH^XY, IH^Xipy (5) 


Here I denotes both charged leptons and neutrinos, and h denotes both neutral and charged 
higgses, and likewise for their superpartners I and H. The above processes are suppressed by M* 
and therefore not strong enough to bring X into equilibrium. Rather, since these interactions 
are extremely feeble, the abundance of X gradually builds up via the process of freeze-in [52] 
as long as the processes remain kinematically feasible. Given the nonrenormalizable operators 
that leads to these interactions, the interaction cross sections scale as ~ s/M^, where s is the 
center of mass energy of the annihilation process, and the production rate is proportional to 
the temperature of the Universe, being the greatest at the earliest times. 

The same processes also result in freeze-in abundances of Y, ipx, and yfy- Assuming 
> my > my and > mx, these particles then decay via 

t[jx ^Ylh, Y ^XIH, ijy^Xlh, (6) 


converting the abundances of U, ipy^ and into X abundance. Taking all these contributions 
into account, we calculate the relic abundance of X to be approximately (in agreement with 
previously derived results in [53H56]) 


~ 0.12 


mx 

10 PeV 



T, 


RH 


1.5 X IQio GeV 


/lO^® GeVy 

V ) 


(7) 


where we have taken a = ai for simplicity. Therefore, with a sufficiently high reheat temper¬ 
ature Thh and appropriate values of a, the PeV scale particle X could compose a significant 
fraction of dark matter. 


2. Decay of X 

Next, we must ensure that X has a lifetime much longer than the age of the Universe and 
the correct decay rate and channels to produce the neutrinos observed at IceGube. We have 
already chosen m^y , my > my, hence the only term in the superpotential Eq.j^that can cause 
X to decay is . Assuming (0)>m<^, the leading decay process is X —>■ coming 
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from the Lagrangian term 


3 


As f (0) 




(8) 


Here we assume, for simplicity, that the decay X —)■ 00 or the decays from mixing with 0 
induced by the corresponding soft term ^X0® that appears after supersymmetry breaking are 
subdominant (however, we have checked that these channels also give similar neutrino spectra 
to that from X —>■ Assuming m^^/mx -C 1, this decay process has a lifetime 

^- 3\2 


Tx 




1.5 X 10“ 
As 


M* 


1016 GeV 


100 PeV 


PeV 

mx 


(9) 


The 00 further decays as 'ijj^ ^ NHu ^ NH^P through an off-shell sterile sneutrino as 
a consequence of the LiHuAfj^ and A/iA/i<h<h terms in the superpotential. The sterile neutrinos 
N then further decay through the standard sterile neutrino decay channels to produce additional 
active neutrinos. 

As the decay lifetime required to £t the IceCube data is r ~ Idol's, Eq.[^ suggests that one 
can obtain the necessary lifetime for reasonable choices of parameters in the model (see section 


the motivation behind the choice of the t/(l)' charge of -|-5 for A. 


Ill A below). Note the role of the 1/M^ suppression in obtaining such a long lifetime; this was 


III. COMPATIBILITY WITH SIGNALS 

In this section we demonstrate the compatibility of the IceCube neutrino and 3.5 keV X-ray 
line signals with the framework described in the previous section. As mentioned in the intro¬ 
duction, one could incorporate neither, one, or both of these into the model with appropriate 
parameter choices. In this section we choose to include both, as a proof of principle that both 
can be incorporated simultaneously into the framework. To demonstrate this, we work with a 
specihc choice of parameters, which are listed in Table |TT] the active and sterile neutrino masses 
and relic abundances of various dark matter components that result from these choices are also 
listed. As stressed in the Introduction, these are not best-£t points resulting from some scan 
but simply a judicious choice of parameters to achieve the desired results. 

The choice tan 0 = 2 is compatible with the measured Higgs mass rrih = 125 GeV with 
PeV scale superpartners. The cutoff scale M* is chosen to be the scale of grand unification 
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Couplings 


/ 

i 


Cij Vi 


3.53 

-2.28 -1.19 X 10“® 

5.82 X 10“® 

a 

= 0.007 

1.02 

-3.54 -1.99 X 10“® 

1.26 

A2 

= 0.0002 

-0.65 

-1.28 3.38 X 10“® ^ 

1.67 




Masses 


mg (eV) 
7.75 X 10-'^ 
0.0087 
0.049 


Dark Matter Properties 


rriNi = 7 V 
= 0.03 (= 25%) 


rus 


7.00 keV 

mx = 7 PeV 

1.50 GeV 

= 2 PeV 

2.00 GeV 

m^o = 800 GeV 

mx = 7 PeV 

mjjo = 800 GeV 

nxh^ = 0.03 (= 25%) 

= 0.06 (= 50%) 


rx = 3 X 10^’^ s 


TABLE II. Our choice of couplings in the superpotential (defined in Eq[^ and the resulting neutrino 
masses and dark matter properties, mg and mg denote the three active and sterile neutrino masses 
respectively. Along with these choices, we have set (</>) = 110 PeV, M* = Mqut (= 10^® GeV), 
tan/3 = 2, and Trh = 10^® GeV. 


Mgut (= 10^® GeV), so the framework is expected to be embedded in a grand unified theory. 
With ( 0 ) = 110 PeV, the specihed values of Qj and rji set the masses of sterile and active 
neutrinos (m^ and respectively) via the see-saw mechanism. It can be seen that the entries 
are mostly 0{1), except for the third column of Qj and the hrst entry in r]i, which are 0(10“®); 


as mentioned in Section IIB , this hierarchy is made inevitable by the need for Ni to be at the 
keV scale and Ai^ 2,3 to be at the GeV scale for consistency with cosmology. It is worth noting 
that although 0(10“®) seems unnaturally small, such a small number is already realized in 
nature in the form of the electron Yukawa coupling. With these choices, the lightest sterile 
neutrino W has a mixing angle of sin^(20) ~ 4 x 10“^® and accounts for ~ 25% of the dark 
matter abundance. 
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FIG. 1. The mass spectrum of our model. The left box displays the heavy PeV scale states as well 
as the TeV scale Higgsinos. The right box shows the active and sterile neutrino mass eigenstates. 
Particles that form some fraction of dark matter are denoted by dashed blue lines. 


As mentioned earlier, the Higgsino is chosen to be the LSP. We take its mass to be 800 GeV, 
so that it makes up about half of the dark matter. The scalar X with mass 7 PeV makes up 
the remaining fraction. Its abundance is controlled by the parameter a (for simplicity, we have 
taken a universal value ctj = a) and the reheat temperature (which we have set to 10^° 
GeV). Likewise, its decay rate is controlled by the parameter A 2 . We see that the correct relic 
density and decay rate can be obtained for fairly reasonable values of these couplings. 

The spectrum of masses of various particles in the model relevant for this paper are plotted 
in Figure[^ with particles that contribute to dark matter highlighted in blue. We have an 
interesting scenario where dark matter is a mix of three components at very different scales: 
a 7 keV sterile neutrino, a stable 800 GeV Higgsino, and a 7 PeV long-lived scalar. For the 
dark matter distribution we adopt the Navarro-Frenk-White (NFW) prohle [57], which was also 
found to be compatible with models of warm [58] and mixed [59] dark matter, as well as the 
3.5 keV line [H]. 

With these parameter choices, we now examine, in turn, the two signals of interest. 
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A. PeV Neutrinos at IceCube 


In this section we study the spectrum of neutrinos at IceCube from the decay of X with the 
parameter choices above. For decaying dark matter, both galactic and extragalactic contribu¬ 
tions are important. The galactic contribution to the neutrino flux d^a/dE^^, where a denotes 
the neutrino flavor, is given by 


dE„ 


= K 


gal 


-13 


1.5 X 10 


cm^ sr s 


1028 ; 

TX 


PeV\ 1 dN^ 


mx J N dEi, 


( 10 ) 


where k = 0.25 is the fraction of dark matter made up by X, l/N{dNa/dEy^) is the normalized 
neutrino energy spectrum from X decay, and tx and mx are, respectively, the lifetime and mass 
of X. Likewise, the extragalactic contribution is [21 EO] 


dE„ 


= K 


ex-gal 


-13 


2.5 X 10 


cm^ sr s 


1028 c 

Tx 


Pev\ r 

mx J Ji 


dy 


dNry 


y 


-3/2 


Nd{E,^y) ^1 + (flA/flM) y-^ ’ 

( 11 ) 


where Qm = 0.315, Ha = 0.685, and y = 1 + z, where 2 ; is the redshift. The expected number 
of neutrino events at IceCube in the energy bin between Ei and E 2 is given by 

d^r 


N{E,,E2) = 


_ /•E 2 




dE„ 


( 12 ) 


^OL / gal -h ex-gal 


where T = 988 days is the total exposure time [T] and A^Q{Ey) is the IceCube effective area for 
neutrino flavor a, taken from |61j . 

We entered the model described in Section into PYTHIA 8.2 [62] and generated the 


neutrino spectrum, dNa/dEu^ from X decay following the decay chain specihed in Section IIC 


The various decay channels and rates for the (P(keV-GeV) mass sterile neutrinos are listed 
in the appendix of [H^. The number of neutrino events in each IceCube bin thus predicted 
following the above procedure is shown in Figure]^ in blue. The IceCube data points (black 
dots), error bars, and background (green) are taken from [1]. The background events come 
from cosmic ray air showers, which produce vr/iC, and expected charmed mesons, which decay 
to muons and neutrinos; the background reported here is the sum of the average background 
from atmospheric muons and neutrinos, plus the 90% conhdence level upper limit for neutrinos 


^ The coefficient in the first parenthesis differs from that in [2] since we use an NFW profile instead of an 
Einasto profile. 
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FIG. 2. The expected number of neutrino events at IceCube from the decay of X (blue), together 
with IceCube data points and error bars (black) and atmospheric background (green) from [T], and 
the total number of signal+background events (red). 

from charmed decays |1]. The total number of events expected (signal+background) for the 
reported exposure is shown in red. 

From the plot, it is clear that the measurements in the lower energy bins are completely 
consistent with atmospheric background, whereas the higher energy bins (200 TeV and above) 
show a clear deviation from what is expected from background only, indicating an excess of 
neutrino events that requires some explanation. This excess has four salient features: 

(i) an excess in two bins at 200 — 400 TeV, 

(ii) no events in the two bins covering 500 TeV - 1 PeV, 

(iii) three PeV scale events, two at 1 PeV and one at 2 PeV, and 

(iv) no events above 2 PeV. 

An additional astrophysical power-law contribution could address (i), but clearly cannot explain 
the three neutrinos with PeV energies (iii) [2]. 

Models of PeV dark matter in the literature claim to be able to explain these features by 
employing a two-body decay channel that includes a neutrino to give (iii) as well as a secondary 
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decay channel that gives softer neutrinos to explain (i) [2HI2]- In contrast, our model does not 
have a direct two-body decay channel into neutrinos, hence there is no sharp feature that peaks 
at PeV energies. Instead, since neutrinos are produced via two- or three- step decay chains, each 
involving multiple decay products, the neutrino spectrum is essentially flat, with a dropoff at 
approximately half the dark matter mass (blue curve in Figurej^. The flat spectrum allows us 
to generate a signal contribution that satisfactorily addresses both (i) and (hi), but at the cost 
of disagreeing with (ii). However, since less than 2 events are predicted in each of the two bins in 
(ii), the disagreement is not too fatal and can be attributed to a possible downward fluctuation 
of the signal. Hence the model predicts that events should appear in these bins when more 
data is collected. Finally, although the DM mass is chosen to have the spectrum drop off after 
the 2 PeV bin, hence explaining (iv), we note that somewhat heavier masses would still predict 
less than 1 event in each of the bins higher than 2 PeV and hence could remain compatible 
with the current data. Consequently, possible future measurements of events at energies higher 
than 2 PeV need not necessarily be incompatible with a dark matter explanation. 

Additional Constraints 

In addition to neutrinos, the dark matter decay products can also give visible signatures in 
other channels. The most important of these is gamma rays. The from DM decays produce 
energetic photons due to inverse Compton scattering and synchrotron radiation. Likewise, since 
the Universe is opaque to gamma rays with energies above a TeV, high energy gamma rays 
produce pairs through interactions with the interstellar radiation held. Such cascades from 
high energy products from DM decay therefore produce a population of gamma rays between 
0{1) GeV and 0(100) GeV [M]- Following [7], we verify the compatibility of the DM decay 
process with gamma-ray bounds from the Fermi-LAT measurement of the isotropic diffuse 
gamma-ray background [65] by considering the integrated energy density measured by the 
Fermi-LAT between Ei ~ 0(1) GeV and E 2 ~ 0(100) GeV [7] 

= — [ E-y—CdE^ ^ 4: A X eV/cm^ . (13) 

c J El d,E^ 

The total energy density in photons and from the decay of X, using the output from 
PYTHIA, is calculated to be ~ 4 x 10“® eV/cm^. This can be interpreted as the maximum 
amount of energy that can be deposited in the diffuse gamma-ray background; since this is 
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well below the energy density measnred by Fermi quoted above, we conclude that the DM 


decay process is not in tension with gamma-ray measurements Likewise, we note that since 
the DM decays exclusively to leptons, it is unlikely to be in tension with the recent Ams- 02 
measurements of the p/p ratio [6^ (from extrapolating Figure 5 of [68], which shows the bound 


on the lifetime of a DM particle decaying to bb from antiproton constraints, to 0(10) PeV, it 


can be seen that even a hadronically decaying particle with a lifetime of ~ lO^^s would be well 


within the antiproton limits). 

B. 3.5 keV X-ray Signal 

Several papers have interpreted this signal as the decay of a ~ 7 keV sterile neutrino dark 


matter component psi o ESI EMS). A keV scale sterile neutrino can be produced through 
the Dodelson-Widrow (DW) mechanism [SU] with an approximate relic abundance 



(14) 


This is a consequence of non-resonant oscillation due to the mixing between the active and 
sterile sectors. However, sterile neutrinos produced through the DW mechanism accounting for 
all of dark matter has now been robustly ruled out for all masses based on constraints from 
X-ray bounds [221 - I29] and Lyman-alpha data [aiiiasiiM] ; see m for a summary. 

This problem can be evaded in several ways. A lepton asymmetry in the early Universe can 
lead to resonant production of sterile neutrinos, resulting in a colder distribution that can evade 
Lyman-alpha bounds; this is known as the Shi-Fuller mechanism [HI]. Another approach is to 
consider sterile neutrino production from the decays of a singlet Higgs boson from an extended 
Higgs sector [85l [HE], which can produce all of dark matter from the freeze-in mechanism [52] 
(for a discussion in the context of the 3.5 keV line, see [70]). 

Likewise, the constraints are avoided if sterile neutrinos make up only a fraction of dark 
matter. While it is trivial to rescale the X-ray constraints to account for a smaller fraction 
of dark matter, a reinterpretation of the Lyman-alpha constraint is not straightforward and 
requires numerical simulations. An analysis in [51], for instance, showed that > 5 keV 
^ See [66j for more stringent (but model-dependent) constraints from analyzing the spectrum rather than the 
integrated energy. 
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warm component constituting < 60% of the total dark matter abundance, where the rest of 
the dark matter is made up of a cold component, is consistent with all constraints ISIli a 
follow-up study |88] found similar results. In a warm plus cold dark matter setup, [59] found 
that the 3.5 keV X-ray signal could be explained with a 7 keV sterile neutrino produced from 
DW that made up 10 — 60% of dark matter. We construct our theory to map on to one of 
the benchmark points in [59], which were shown to satisfy the relevant constraints. In [591, the 
signal was found to be compatible with a non-resonantly produced 7 keV sterile neutrino with 
mixing angle sin^(26*) ~ 4 x 10“^°, making up ~ 25% of dark matter; we reproduce these values 
with our choice of parameters in Table [IT] thereby incorporating the 3.5 keV X-ray line in our 
framework. 

A mixture of warm and cold components for dark matter offers several advantages. Ref. 
[59] found that this scenario was still compatible with an NFW prohle and could resolve the 
missing satellite problem [891199] . In addition, we note that the larger free-streaming length of 
the non-resonantly produced sterile neutrino, which is a warm dark matter component, could 
result in it being underabundant in dwarf galaxies (which would be made up mostly of the 
cold component), possibly providing an explanation of the non-observation of the signal in the 
stacked analysis of dwarf spheroidals presented in [TH] ; establishing this would, however, require 
detailed numerical simulations that are beyond the scope of this paper. 


IV. SUMMARY 

In this paper we have attempted to incorporate two recent potential hints of new physics, the 
PeV neutrinos at IceCube and the 3.5 keV X-ray line, into a broader, independently motivated 
framework of a PeV scale supersymmetric neutrino sector. 

• The right handed sterile neutrinos are expected to be charged under some new symmetry 
(e.g., 77(1)') beyond the SM gauge group, which enables light masses at phenomenologi¬ 
cally interesting scales (keV-GeV). In order to successfully realize the desired active and 
sterile neutrino masses without unnaturally small parameters, this symmetry must be 
broken by a PeV scale vev, which corresponds to a desired scale in some approaches to 
supersymmetry breaking. 
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• The lightest sterile neutrino can have a mass of 7 keV and be produced non-resonantly 
through the DW mechanism to form a fraction of dark matter; its decay can explain the 
3.5 keV X-ray line. 

• The neutrino sector can be extended to include additional helds with similar operators 
to the ones that give rise to neutrino masses. This can give a PeV scale dark matter 
component whose relic abundance is set by freeze-in processes, and an appropriate charge 
under the U{1)' symmetry makes it long-lived and gives a decay spectrum into neutrinos 
consistent with the high energy events observed at IceCube. 

Although we have worked in a setup where dark matter consists of 25% sterile neutrinos, 
25% PeV scalar X, and 50% Higgsino LSP, it should be clear that the dark matter composition 
can be completely different from this admixture, given that the PeV supersymmetric sector is 
likely far more complicated, and neither of the two signals are therefore necessary ingredients 
of the extended neutrino sector we have discussed. Since what is presented here only serves as 
an illustrative proof of principle and can incorporate several variations, we did not attempt to 
perform detailed hts or delve into model-building details beyond what was necessary, although 
such exercises might be warranted once more data is collected in the future, and can shed 
light on the underlying model. These should not distract from the most central ideas of the 
paper that these two signals, at such different energy scales, £t rather naturally into a broader 
particle physics framework and could be the hrst indications of a rich, hitherto unexplored 
neutrino sector. 

Acknowledgements: The authors are supported in part by the DoE under grants DE- 
SC0007859 and DE-SC0011719. BS also thanks the Center for Theoretical Underground Physics 
and Related Areas (CETUP*), where part of this work was completed, for its hospitality. 


[1] M. Aartsen et al. (IceCube), Phys.Rev.Lett. 113, 101101 (2014), arXiv: 1405.5303 [astro-ph.HE], 

[2] C. S. Fong, H. Minakata, B. Panes, and R. Z. Funchal, JHFP 1502, 189 (2015), arXiv:1411.5318 
[hep-ph], 

[3] Y. Fma, R. Jinno, and T. Moroi, Phys.Lett. B733, 120 (2014), arXiv:1312.3501 [hep-ph], 


17 




[4] Y. Ema, R. Jinno, and T. Moroi, JHEP 1410, 150 (2014), arXiv:1408.1745 [hep-ph], 

[5] B. Eeldstein, A. Kusenko, S. Matsumoto, and T. T. Yanagida, Phys.Rev. D88, 015004 (2013), 
arXiv:1303.7320 [hep-ph], 

[6] A. Esmaili and P. D. Serpico, JCAP 1311, 054 (2013), arXiv: 1308.1105 [hep-ph], 

[7] A. Esmaili, S. K. Kang, and P. D. Serpico, JCAP 1412, 054 (2014), arXiv:1410.5979 [hep-ph], 

[8] Y. Bai, R. Lu, and J. Salvado, (2013), arXiv:1311.5864 [hep-ph], 

[9] A. Bhattacharya, M. H. Reno, and I. Sarcevic, JHEP 1406, 110 (2014), arXiv:1403.1862 [hep-ph], 

[10] T. Higaki, R. Kitano, and R. Sato, JHEP 1407, 044 (2014), arXiv:1405.0013 [hep-ph], 

[11] A. Bhattacharya, R. Gandhi, and A. Gupta, JCAP 1503, 027 (2015), arXiv:1407.3280 [hep-ph], 

[12] C. Rott, K. Kohri, and S. C. Park, (2014), arXiv: 1408.4575 [hep-ph], 

[13] E. Bulbul, M. Markevitch, A. Foster, R. K. Smith, M. Loewenstein, et al., Astrophys.J. 789, 13 
(2014), arXiv: 1402.2301 [astro-ph.CO]. 

[14] A. Boyarsky, O. Ruchayskiy, D. lakubovskyi, and J. Franse, Phys.Rev.Lett. 113, 251301 (2014), 
arXiv:1402.4119 [astro-ph.CO], 

[15] T. E. Jeltema and S. Profumo, (2014), arXiv:1408.1699 [astro-ph.HE], 

[16] A. Boyarsky, J. Franse, D. lakubovskyi, and O. Ruchayskiy, (2014), arXiv: 1408.4388 [astro- 
ph.CO], 

[17] E. Bulbul, M. Markevitch, A. R. Foster, R. K. Smith, M. Loewenstein, et al, (2014), 
arXiv: 1409.4143 [astro-ph.HE], 

[18] D. Malyshev, A. Neronov, and D. Eckert, Phys.Rev. D90, 103506 (2014), arXiv: 1408.3531 [astro- 
ph.HE] , 

[19] M. E. Anderson, E. Churazov, and J. N. Bregman, (2014), arXiv:1408.4115 [astro-ph.HE], 

[20] S. Riemer-Sorensen, (2014), arXiv:1405.7943 [astro-ph.CO], 

[21] A. Boyarsky, J. Franse, D. lakubovskyi, and O. Ruchayskiy, (2014), arXiv: 1408.2503 [astro- 
ph.CO], 

[22] A. Boyarsky, A. Neronov, O. Ruchayskiy, M. Shaposhnikov, and 1. Tkachev, Phys.Rev.Lett. 97, 
261302 (2006), arXiv:astro-ph/0603660 [astro-ph], 

[23] A. Boyarsky, J. Nevalainen, and O. Ruchayskiy, Astron.Astrophys. 471, 51 (2007), arXiv:astro- 
ph/0610961 [astro-ph]. 


18 












[24] C. R. Watson, J. F. Beacom, H. Yuksel, and T. P. Walker, Phys.Rev. D74, 033009 (2006), 
arXiv:astro-ph/0605424 [astro-ph]. 

[25] A. Boyarsky, A. Neronov, O. Ruchayskiy, and M. Shaposhnikov, Mon.Not.Roy.Astron.Soc. 370, 
213 (2006), arXiv:astro-ph/0512509 [astro-ph], 

[26] A. Boyarsky, D. lakubovskyi, O. Ruchayskiy, and V. Savchenko, Mon.Not.Roy.Astron.Soc. 387, 
1361 (2008), arXiv:0709.2301 [astro-ph], 

[27] A. Boyarsky, D. Malyshev, A. Neronov, and O. Ruchayskiy, Mon.Not.Roy.Astron.Soc. 387, 1345 
(2008), arXiv:0710.4922 [astro-ph], 

[28] C. R. Watson, Z.-Y. Li, and N. K. Policy, JCAP 1203, 018 (2012), arXiv:1111.4217 [astro-ph.CO], 

[29] S. Horiuchi, P. J. Humphrey, J. Onorbe, K. N. Abazajian, M. Kaplinghat, et al., Phys.Rev. D89, 
025017 (2014), arXiv:1311.0282 [astro-ph.CO], 

[30] K. Abazajian, Phys.Rev. D73, 063513 (2006), arXiv:astro-ph/0512631 [astro-ph], 

[31] U. Seljak, A. Makarov, P. McDonald, andH. Trac, Phys.Rev.Lett. 97, 191303 (2006), arXiv:astro- 
ph/0602430 [astro-ph], 

[32] T. Asaka, M. Laine, and M. Shaposhnikov, JHEP 0701, 091 (2007), arXiv:hep-ph/0612182 
[hep-ph], 

[33] M. Viel, J. Lesgourgues, M. G. Haehnelt, S. Matarrese, and A. Riotto, Phys.Rev.Lett. 97, 071301 
(2006), arXiv:astro-ph/0605706 [astro-ph], 

[34] A. Boyarsky, J. Lesgourgues, O. Ruchayskiy, and M. Viel, JCAP 0905, 012 (2009), 
arXiv:0812.0010 [astro-ph], 

[35] E. Polisensky and M. Ricotti, Phys.Rev. D83, 043506 (2011), arXiv: 1004.1459 [astro-ph.CO], 

[36] M. Viel, G. D. Becker, J. S. Bolton, and M. G. Haehnelt, Phys.Rev. D88, 043502 (2013), 
arXiv:1306.2314 [astro-ph.CO], 

[37] S. B. Roland, B. Shakya, and J. D. Wells, (2014), arXiv:1412.4791 [hep-ph], 

[38] J. D. Wells, (2003), arXiv:hep-ph/0306127 [hep-ph], 

[39] N. Arkani-Hamed and S. Dimopoulos, JHEP 0506, 073 (2005), arXiv:hep-th/0405159 [hep-th], 

[40] G. Giudice and A. Romanino, Nucl.Phys. B699, 65 (2004), arXiv:hep-ph/0406088 [hep-ph], 

[41] J. D. Wells, Phys.Rev. D71, 015013 (2005), arXiv:hep-ph/0411041 [hep-ph], 

[42] G. F. Giudice and A. Strumia, Nucl.Phys. B858, 63 (2012), arXiv: 1108.6077 [hep-ph]. 


19 
















[43] N. Arkani-Hamed, A. Gupta, D. E. Kaplan, N. Weiner, and T. Zorawski, (2012), arXiv: 1212.6971 
[hep-ph], 

[44] A. Arvanitaki, N. Craig, S. Dimopoulos, and G. Villadoro, JHEP 1302, 126 (2013), 
arXiv:1210.0555 [hep-ph], 

[45] M. R. Lovell, V. Eke, C. S. Erenk, L. Gao, A. Jenkins, et al., Mon.Not.Roy.Astron.Soc. 420, 2318 
(2012), arXiv: 1104.2929 [astro-ph.CO]. 

[46] M. Boylan-Kolchin, J. S. Bullock, and M. Kaplinghat, Mon.Not.Roy.Astron.Soc. 422, 1203 
(2012), arXiv:111 1.2048 [astro-ph.CO]. 

[47] A. Kusenko, Phys.Rept. 481, 1 (2009), arXiv:0906.2968 [hep-ph], 

[48] P. Hernandez, M. Kekic, and J. Lopez-Pavon, Phys.Rev. D90, 065033 (2014), arXiv: 1406.2961 
[hep-ph], 

[49] A. C. Vincent, E. F. Martinez, P. Hernandez, M. Lattanzi, and O. Mena, (2014), arXiv: 1408.1956 
[astro-ph.CO]. 

[50] S. Dodelson and L. M. Widrow, Phys.Rev.Lett. 72, 17 (1994), arXiv:hep-ph/9303287 [hep-ph], 

[51] K. Griest and M. Kamionkowski, Phys.Rev.Lett. 64, 615 (1990). 

[52] L. J. Hall, K. Jedamzik, J. March-Russell, and S. M. West, JHEP 1003, 080 (2010), 
arXiv:0911.1120 [hep-ph], 

[53] F. Elahi, C. Kolda, and J. Unwin, (2014), arXiv:1410.6157 [hep-ph], 

[54] A. Kusenko, F. Takahashi, and T. T. Yanagida, Phys.Lett. B693, 144 (2010), arXiv:1006.1731 
[hep-ph], 

[55] M. Blennow, E. Fernandez-Martinez, and B. Zaldivar, JCAP 1401, 003 (2014), arXiv:1309.7348 
[hep-ph], 

[56] S. Khalil and O. Seto, JCAP 0810, 024 (2008), arXiv:0804.0336 [hep-ph], 

[57] J. F. Navarro, C. S. Erenk, and S. D. White, Astrophys.J. 462, 563 (1996), arXiv:astro- 
ph/9508025 [astro-ph], 

[58] M. R. Lovell, C. S. Erenk, V. R. Eke, A. Jenkins, L. Gao, et al., Mon.Not.Roy.Astron.Soc. 439, 
300 (2014), arXiv: 1308.1399 [astro-ph.CO], 

[59] A. Harada, A. Kamada, and N. Yoshida, (2014), arXiv: 1412.1592 [astro-ph.CO], 

[60] A. Ibarra and D. Tran, Phys.Rev.Lett. 100, 061301 (2008), arXiv:0709.4593 [astro-ph], 


20 
















[61] M. Aartsen et al. (IceCube Collaboration), Phys.Rev.Lett. Ill, 021103 (2013), arXiv: 1304.5356 
[astro-ph.HE], 

[62] T. Sjstrand, S. Ask, J. R. Christiansen, R. Corke, N. Desai, et al., Comput.Phys.Commun. 191, 
159 (2015), arXiv;1410.3012 [hep-ph], 

[63] D. Gorbunov and M. Shaposhnikov, JHEP 0710, 015 (2007), arXiv:0705.1729 [hep-ph], 

[64] D. V. Semikoz and G. Sigl, JCAP 0404, 003 (2004), arXiv:hep-ph/0309328 [hep-ph], 

[65] M. Ackermann et al. (Fermi-LAT), Astrophys.J. 799, 86 (2015), arXiv:1410.3696 [astro-ph.HE], 

[66] K. Murase, R. Laha, S. Ando, and M. Ahlers, (2015), arXiv:1503.04663 [hep-ph], 

[67] “Ams-02 collaboration, talks at the ‘AMS Days at CERN’, 15-17 april 2015.”. 

[68] G. Giesen, M. Boudaud, Y. Genolini, V. Poulin, M. Cirelli, et al, (2015), arXiv: 1504.04276 
[astro-ph.HE], 

[69] K. N. Abazajian, Phys.Rev.Lett. 112, 161303 (2014), arXiv: 1403.0954 [astro-ph.CO]. 

[70] A. Merle and A. Schneider, (2014), arXiv: 1409.6311 [hep-ph], 

[71] N. Haba, H. Ishida, and R. Takahashi, Phys.Lett. B743, 35 (2015), arXiv:1407.6827 [hep-ph], 

[72] W. Rodejohann and H. Zhang, Phys.Lett. B737, 81 (2014), arXiv: 1407.2739 [hep-ph], 

[73] A. Abada, G. Arcadi, and M. Lucente, (2014), 10.1088/1475-7516/2014/10/001, arXiv: 1406.6556 
[hep-ph], 

[74] H. Ishida and H. Okada, (2014), arXiv: 1406.5808 [hep-ph], 

[75] D. J. Robinson and Y. Tsai, Phys.Rev. D90, 045030 (2014), arXiv:1404.7118 [hep-ph], 

[76] S. Chakraborty, D. K. Ghosh, and S. Roy, JHEP 1410, 146 (2014), arXiv: 1405.6967 [hep-ph], 

[77] A. Adulpravitchai and M. A. Schmidt, JHEP 1501, 006 (2015), arXiv: 1409.4330 [hep-ph], 

[78] S. Patra and P. Pritimita, (2014), arXiv: 1409.3656 [hep-ph], 

[79] M. Frigerio and C. E. Yaguna, (2014), arXiv: 1409.0659 [hep-ph], 

[80] T. Tsuyuki, Phys.Rev. D90, 013007 (2014), arXiv:1403.5053 [hep-ph], 

[81] S. K. Kang and A. Patra, (2014), arXiv:1412.4899 [hep-ph], 

[82] H. Ishida, K. S. Jeong, and F. Takahashi, Phys.Lett. B732, 196 (2014), arXiv:1402.5837 [hep-ph], 

[83] F. Bezrukov and D. Gorbunov, Phys.Lett. B736, 494 (2014), arXiv: 1403.4638 [hep-ph], 

[84] X.-D. Shi and G. M. Fuller, Phys.Rev.Lett. 82, 2832 (1999), arXiv:astro-ph/9810076 [astro-ph], 

[85] A. Kusenko, Phys.Rev.Lett. 97, 241301 (2006), arXiv:hep-ph/0609081 [hep-ph]. 


21 
















[86] K. Petraki and A. Kusenko, Phys.Rev. D77, 065014 (2008), arXiv:0711.4646 [hep-ph], 

[87] A. Boyarsky, O. Ruchayskiy, and M. Shaposhnikov, Ann.Rev.Nucl.Part.Sci. 59, 191 (2009) 
arXiv:0901.0011 [hep-ph], 

[88] A. Boyarsky, J. Lesgourgues, O. Ruchayskiy, and M. Viel, Phys.Rev.Lett. 102, 201304 (2009) 
arXiv:0812.3256 [hep-ph], 

[89] B. Moore, S. Ghigna, F. Governato, G. Lake, T. R. Quinn, et al, Astrophys.J. 524, L19 (1999) 
arXiv:astro-ph/9907411 [astro-ph]. 

[90] A. V. Kravtsov, Adv.Astron. 2010, 281913 (2010), arXiv:0906.3295 [astro-ph. CO], 


22 



